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Abstract. This paper proposes an optimal allocation problem with ramified transport 
technology in a spatial economy. Ramified transportation is used to model the transport 
economy of scale in group transportation observed widely in both nature and efficiently 
designed transport systems of branching structures. The ramified allocation problem aims 
at finding an optimal allocation plan as well as an associated optimal allocation path 
to minimize overall cost of transporting commodity from factories to households. This 
problem differentiates itself from existing ramified transportation literature in that the 
distribution of production among factories is not fixed but endogenously determined as 
observed in many allocation practices. It's shown that due to the transport economy 
of scale in ramified transportation, each optimal allocation plan corresponds equivalently 
to an optimal assignment map from households to factories. This optimal assignment 
map provides a natural partition of both households and allocation paths. We develop 
methods of marginal transportation analysis and projectional analysis to study properties 
of optimal assignment maps. These properties are then related to the search for an optimal 
assignment map in the context of state matrix. 



1. Introduction 



One of the lasting interests in economics is to study optimal resource allocation in a 
spatial economy. For instance, the well known Monge-Kantorovich transport problem aims 
at finding an efficient allocation plan or map for transporting some commodity from fac- 
tories to households. This problem was pioneered by Monge 19 and advanced fully by 



Kantorovich [13 1| who won the Nobel prize in economics in 1975 for his seminal work on 
optimal allocation of resources. Recent advancement of this problem in mathematics can be 



found in Villani 22 , 23 and references therein. Monge-Kantorovich problem has also been 
applied to study other related economic problems, e.g., spatial firm pricing (Buttazzo and 
Carlier [4j), principal-agent problem (Figalli, Kim and McCann |10|), hedonic equilibrium 
models (Chiappori, McCann and Nesheim |7j; Ekeland |9|), matching and partition in labor 
market (Carlier and Ekeland [5]; McCann and Trokhimtchouk [l8]). 

Recently, a new research field known as Ramified Optimal Transportation has grown 
out of Monge-Kantorovich problem. Representative studies can be found for instance 
in Gilbert 12 , Xia [24-30 , Maddalena, Solimini and Morel 16 , Bernot, Caselles and 



Morel [1|[2], Brancolini, Buttazzo and Santambrogio |3|, Santambrogio |21| , Devillanova 
and Solimini [8] , Xia and Vershynina [3j] . Ramified optimal transport problem studies how 
to find an optimal transport path from sources to targets as shown in Figure [TJ Different 
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from the standard Monge-Kantorovich problem where the transportation cost is solely de- 
termined by a transport plan or map, the transportation cost in the ramified transport 
problem is determined by the actual transport path which transports the commodity from 
sources to targets. Ramified transportation indeed formally formulates the concept of trans- 
port economy of scale in group transportation observed widely in both nature (e.g. trees, 
blood vessels, river channel networks, lightning) and efficiently designed transport systems 
of branching structures (e.g. railway configurations and postage delivery networks). An ap- 
plication of ramified optimal transportation in economics can be found in Xia and Xu [32] , 
which showed that a well designed ramified transport system can improve the welfare of 
consumers in the system. 

In this paper, we propose an optimal resource allocation problem where a planner chooses 
both an optimal allocation plan as well as an associated optimal transport path using ram- 
ified transport technology. In both Monge-Kantorovich and ramified transport problems, 
one typically assumes an exogenous fixed distribution in both sources and targets. How- 
ever, in many resource allocation practices, the distribution in either sources or targets is 
not pre-determined but rather determined endogenously. For instance, in a production al- 
location problem, suppose there are k factories and i households located in different places 
in some area. The demand for some commodity from each household is fixed. Nevertheless, 
the allocation of production among factories is not pre-determined but rather depends on 
the distribution of demand among households as well as their relative locations to factories. 
A planner needs to make an efficient allocation plan of production over these k factories 
to meet given demands from I households. Under ramified optimal transportation, the 
transportation cost of each production plan is determined by an associated optimal trans- 
port path from factories to households. Consequently, the planner needs to find an optimal 
production plan as well as an associated optimal transport path to minimize overall cost 
of distributing commodity from factories to households. Another example of similar nature 
exists in the following storage arrangement problem. Suppose there are k warehouses and £ 
factories located in different places in some area. Each factory has already produced some 
amount of commodity. However, the assignment of commodity among warehouses is not 
pre-determined but instead relies on the distribution of production among factories as well 
as relative locations between factories and warehouses. Similarly, a planner needs to make 
an efficient storage arrangement as well as an associated optimal transport path for storing 
the produced commodity in the given k warehouses with minimal transportation cost. 

Problem of this category is formulated as the ramified optimal allocation problem in 
Section 2. Throughout the following context, we will focus our discussion on the scenario of 
the production allocation problem. Little additional effort is needed to interpret results for 
other scenarios. We start with modeling a transport path from factories to households as a 
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weighted directed graph, where the transportation cost on each edge of the graph depends 
linearly on the length of the edge but concavely on the amount of commodity moved on 
the edge. The motivation of concavity of the cost functional on quantity comes from the 
observation of transport economy of scale in group transportation. The more concave is 
the cost functional or the greater is the magnitude of transport economy of scale, the more 
efficient is to transport commodity in larger groups. We define the cost of an allocation 
plan as the minimum transportation cost of a transport path compatible with this plan. A 
planner needs to find an efficient allocation plan such that demands from households will be 
met in a least cost way. In this problem, the distribution of production over factories is not 
pre-determined as in Monge-Kantorovich or ramified transport problems, but endogenously 
determined by the distribution of demands from households as well as their relative locations 
to factories. 

We prove the existence of the ramified allocation problem in Section 3. It's shown that 
due to the transport economy of scale in ramified transportation, under any optimal allo- 
cation plan, no two factories will be connected on any associated optimal allocation path. 
Consequently, any optimal allocation path can be decomposed into a set of mutually dis- 
joint transport paths originating from each factory. As a result, each household will receive 
her commodity from only one factory under any optimal allocation plan. It implies that 
each optimal allocation plan corresponds equivalently to an optimal assignment map from 
households to factories. Thus, solving the ramified optimal allocation problem is equivalent 
to finding an optimal assignment map. This optimal assignment map is shown to provide 
a partition not only in households but also in the associated allocation path according to 
the factories. 

Because of the equivalence between the optimal allocation plan and assignment map, 
we can instead focus attention on studying the properties of optimal assignment maps in 
the ramified optimal allocation problem. In Section 4, we develop a method of marginal 
transportation analysis to study properties of optimal assignment maps. This method 
extends the standard marginal analysis in economics into the analysis for transport paths. 
It builds upon an intuitive idea that a marginal change on an optimal allocation path 
should not reduce the existing minimal transportation cost. Using this method, we develop 
a criterion which relates the optimal assignment of a household with her relative location 
to factories and other households, as well as her demand and productions at factories. In 
particular, it is shown that each factory has a nearby region such that a household living at 
this region will be assigned to the factory, where the size of this region depends positively 
on the demand of the household. In this case, the planner takes advantage of relative 
spatial locations between households and factories. Also, if an optimal assignment map 
assigns a household to some factory, then this household has a neighborhood area such that 
any household with a smaller demand living in this area will also be assigned to the same 
factory. Here, the planner utilizes the benefit in group transportation due to transport 
economy of scale embedded in ramified transportation. The role of spatial location and 
group transportation in resource allocation is further studied in Section 5 by a method of 
projectional analysis. We show that under an optimal assignment map, a household will 
be assigned to some factory only when either she lives close to the factory or she has some 
nearby neighbors assigned to the factory. In particular, there is an "autarky" situation 
when households and factories are located on two disjoint areas lying distant away from 
each other, the demand of households will solely be satisfied from local factories. 
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An important application of the properties of optimal assignment maps is that they can 
shed light on the search for those maps. In Section 6, we develop a search method utilizing 
these properties in a notion of state matrix. A state matrix represents the information set 
of a planner during the search process for an optimal assignment map. Any zero entry u s h 
in the matrix reflects that the planner has excluded the possibility of assigning household 
h to factory s under this map. When a state matrix has exactly one non-zero entry in each 
column, it completely determines an optimal assignment map by those non-zero entries. 
Our search method uses properties about optimal assignment maps to update some non- 
zero entries with zeros in a state matrix. This method is motivated by the observation that 
via group transportation under ramified transport technology, assignment of each household 
has a global effect on the allocation path as well as the associated assignment map. Thus, 
the planner can deduce more information about the optimal assignment map by exploiting 
the existing information embedded in zero entries of a state matrix. Each updated state 
matrix contains more zeros and thus more information than its pre-updated counterpart. 
This method is useful in the search for optimal assignment maps as each updating step 
increases the number of zero entries which in turn reduces the size of the restriction set 
of assignment maps in a large magnitude. In some non-trivial cases, it's shown that this 
method can exactly find an optimal assignment map as desired. 

2. Ramified Optimal Allocation Problem 

In this section, we describe the setting of the optimal allocation model with ramified 
optimal transportation. 

2.1. Ramified Optimal Transportation. In a spatial economy, there are k factories and 
t households located at x = {x\, X2,- Xk} and {2/1,2/2, ■ ■ • , He} in some area X, where A is a 
compact convex subset of a Euclidean space R m . In this model economy, there is only one 
commodity, and each household j = 1, ■ ■ ■ , i has a fixed demand rij > for the commodity. 

For analytical convenience, we first represent households and factories as atomic Radon 
measures. Recall that a Radon measure c on X is atomic if c is a finite sum of Dirac 
measures with positive multiplicities, i.e., 

s 

i=i 

for some integer s > 1 and some points z% G X with a > for each i = 1, ■ ■ ■ , s. The mass 
of c is denoted by 

s 

m(c) :=J2 c i- 

i=i 

We can thus represent the £ households as an atomic measure on X by 

t 

I' Y, "A • (2-1) 

3=1 

For each i = 1, ■ ■ ■ , k, denote mi as the units of the commodity produced at factory i located 
at xi. Then, the k factories can be represented by another atomic measure on X by 

k 

& = ^2mi5 Xt . (2.2) 
i=i 
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In the study of transport problems, we usually assume m (a) = m (b), i.e., 

k I 

i=i j=i 

which simply means that supply equals demand in aggregate. 

Next, we introduce the concept of transport path from a to b as in Xia [24]. 

Definition 2.1. Suppose a andh are two atomic measures on X of equal mass. A transport 
path from a to h is a weighted directed graph G consisting of a vertex set V (G), a directed 
edge set E (G) and a weight function w : E (G) — > (0, +oo) such that {x±, X2, U 
{yi,U2, ■■■,yi} ^ V{G) and for any vertex v £ V(G), there is a balance equation 

{nii, if v = Xi for some i = 1, k 
-rij, if v = yj for some j = !,...,£ (2.3) 
0, otherwise 

where each edge e £ E (G) is a line segment from the starting endpoint e~ to the ending 
endpoint e + . Denote Path (a, b) as the space of all transport paths from a to b. 



Note that the balance equation (2.3) simply means the conservation of mass at each 



vertex. Viewing G as an one dimensional polyhedral chain, equation (2.3) may simply be 
expressed as dG = b — a. 

Now, we consider the transportation cost of a transport path. As observed in both 
nature and efficiently designed transport networks, there exists a transport economy of scale 
underlying group transportation. For this consideration, ramified optimal transport theory 
uses a cost functional depending concavely on quantity and defines the transportation cost 
of a transport path as follows. 

Definition 2.2. For each transport path G G Path (a, b) and any a G [0, 1], the M Q cost 
of G is defined by 

M a (G):= ^ [w (e)f length (e) . (2.4) 

edE{G) 

The parameter a represents the magnitude of transport economy of scale. The smaller 
the a, the more efficient is to move commodity in groups. Ramified optimal transport 
problem studies how to find a transport path to minimize the M a cost, i.e., 

min M a (G) , (2.5) 

GePath(a,b) 

whose minimizer is called an optimal transport path from a to b. An optimal transport 



path has many nice properties. For instance, it contains no cycles by Xia 24, Proposition 
2.1]. Thus, without loss of generality, we assume that all transport paths considered in the 
following context contain no cycles. When a < 1, an optimal transport path is generally of 
branching structure. In the scenario with two sources and one target, a "Y-shaped" path 
is usually more preferable than a "V-shaped" path. 

The following example illustrates the effect of transport economy of scale on optimal 
transport path in a spatial economy with one factory a = 5o located at origin and fifty 
households b = Y^,j=i BU^j °^ ec l ua l demand rij = g^. The locations of these fifty households 
are randomly selected. As seen in Figure [2j when a = 1, the optimal transport path 
is "linear" in the sense that the factory will ship commodity directly to each household. 
When a < 1, the transport path becomes "ramified" as the planner would like commodity 
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Figure 2. Examples of optimal transport paths. 



to be transported in groups in order to utilize the benefit of transport economy of scale. 
Furthermore, by comparing for instance the width of the transport paths for a = 0.75 
and a = 0.25, we observe that the smaller the a, the more likely the commodity will be 
transported in a large scale. 

For any atomic measures a and b on X of equal mass, define the minimum transportation 
cost as 

d a (a, b) := min {M a (G) : G G Path (a, b)} . (2.6) 

As shown in Xia |24|, d a is indeed a metric on the space of atomic measures of equal mass. 
Also, for each A > 0, it holds that 

d a (Aa, Ab) = X a d a (a, b) . 

Without loss of generality, we normalize both a and b to be a probability measure on X, 
i.e., 

k i 

Y j m i = Y j n j = l. (2.7) 

i=l j=l 

2.2. Compatibility between Transport Plan and Path. For the allocation problem, 
a key decision a planner needs to make is about the transport plan from factories to house- 
holds. 



Definition 2.3. Suppose a and b are two atomic probability measures on X as in (2.2), 
(2.1) and (2.1). A transport plan from a to b is an atomic probability measure 



1 = ^2^2 < lij S (x l ,y ] ) 
i=l j=l 

on the product space X x X such that for each i and j, qij > 0, 

k I 
y^qij = Tij and S^ij = mi. 
i=l j=l 

Denote Plan (a, b) as the space of all transport plans from a to b. 



(2.1 



(2.9) 



In a transport plan q, the number qij denotes the amount of commodity received by 
household j from factory i. 

Now, as in Section 7.1 of Xia [24], we want to consider the compatibility between a 
transport path and a transport plan. Let G be a given transport path in Path (a, b). Since 
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G contains no cycles, for each Xi and yj, there exists at most one directed polyhedral curve 
gij on G from x\ to yj. In other words, there exists a list of distinct vertices 

v (9ij) ■= {vh,Vi 2 ,--- ,Vi h } (2.10) 



in V (G) with Xi 



V;, 



Vj 



i h , and each [«i t ,«i t+1 ] is a directed edge in E (G) for each 



t = 1,2, ■ ■ ■ , h — 1. For some pairs of (i,j), such a curve ^ from Xi to j/j may not exist, in 
which case we set = to denote the empty directed polyhedral curve. By doing so, we 
construct a matrix 

9 = (9ij) kx t (2-11) 
with each element of g being a polyhedral curve. For any transport path G G Path (a, b), 
such a matrix g = (g^) is uniquely determined. 

Definition 2.4. Let G £ Path (a, b) be a transport path and q G Plan (a, b) be a transport 
plan. The pair (G, q) is compatible if qij = whenever gij = and 

G = q-g. (2.12) 



Here, equation (2.12) means that as polyhedral chains, 

k i 

G = Y J 12 qi r9ij, 

i=l j=l 

where the product qij ■ gtj denotes that an amount qij of commodity is moved along the 
polyhedral curve g^ from factory i to household j. 






(A) Gi (B) G 2 

Figure 3. Compatibility between transport plan and transport path. 

Roughly speaking, the compatibility conditions check whether a transport plan is real- 
izable by a transport path. Given a transport plan, the planner must design a transport 
path which can support this plan. To see the concept more precisely, let 

13 5 3 

a =-5 xl + -5 X2 and b =-5 yi + -S y2 , 



and consider a transport plan 



g*(«i*i) + g*(iei,«a) + 2^2,2/1) + 4<Wj/ 2 ) g Plan (a,b) . 



(2.13) 
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It is straight forward to see from Figure [3] that q is compatible with G\ but not Gi- This is 
because there is no directed curve g\2 from factory 1 to household 2 in C?2. 

2.3. Ramified Allocation Problem. In the standard transport problems, e.g. Monge- 
Kantorovich or ramified optimal transportation, one typically assumes an exogenous fixed 
distribution of production among factories. In this paper, we consider a scenario where this 
distribution is not fixed but endogenously determined. In other words, the atomic measure 
a which represents the k factories can have varying production level n%i at each factory 
i. This consideration is motivated by our observation that in many allocation practices as 
discussed in the introduction, the distribution of production among factories is not pre- 
determined but rather depends on the distribution of demand among households as well as 
their relative locations to factories. 

Definition 2.5. Let x = {x\,X2,— be a finite subset of X , and b be the atomic proba- 
bility measure representing households defined in (2.1). An allocation plan from x to b is a 
probability measure 

k i 
i=l j=l 

on X x X such that q-ij > for each i, j and 

k 

y j Qij = rij for each j = !,■■■ 
i=l 

Denote Plan [x, b] as the set of all allocation plans from x to b. 

Note that any allocation plan q £ Plan [x, b] corresponds to a transport plan q from a (g) 
to b, where a (q) is the probability measure representing k factories defined as 

k i 
a (^) := ^Z mi *xj, with rxij (q) = y~]gij, i = l,...,k. (2.14) 

i=l j=\ 

In other words, Plan [x, b] is the union of Plan (a, b) among all atomic probability measures 
a supported on x. 

Example 2.1. Any function S : {1, • • • , £} — > {1, • • • , k} determines an allocation plan in 
Plan [x, b] as 

k I 

?S = EE 1iA*iW) With »i = 
i=l j=l 

That is, 



nj, ifi = S(j) 
0, else 



»=E^(«(s W )*)- (2 - 15) 

i=i 

For a given allocation plan, we define the associated transportation cost as follows. 

Definition 2.6. For any allocation plan q € Plan [x, b] and a £ [0, 1), the ramified trans- 
portation cost of q is 

T a (q) := min {M a (G) : G G Path (a (q) , b) , (G,q) compatible}, (2.16) 
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where M a (•) is defined in (2.4\) - An allocation plan q* E Plan [x, b] is optimal if 

T a (q*) < T a (q) for any q G Plan [x, b] . 

For each allocation plan q, as in Xia |24l Proposition 7.3], there exists a path G q G 
Path (a (g) , b) such that G q is compatible with q and 

T Q (g) = M a (G q ) . (2.17) 



Thus, the minimum value in (2.16) is achieved by G q for each q. Now, we are ready to 
define the major problem in this paper. 

Problem 2.1. (Ramified Optimal Allocation Problem) Let X be a compact convex domain 
in M m with the standard norm ||-||. Given a finite subset x = {xi,X2 t — Xk} in X, an atomic 
probability measure h on X defined in (2.1), and a parameter a G [0, 1). Find a minimizer 
ofT a (g) among all allocation plans q G Plan [x, b], i.e., 



min {T Q (q) : q G Plan [x, b]} . 
3. Characterizing Optimal Allocation Plans 



(2.18) 



In this section, we first establish the existence result of the ramified optimal allocation 
problem. It's then shown that any optimal allocation plan corresponds to an optimal 
assignment map from households to factories, which provides a partition in both households 
and transport paths. 



The following proposition proves the existence of Problem 2.1 



Proposition 3.1. The ramified optimal allocation problem (2. IS) has a solution. Moreover, 

:mGifj, (3.1) 

where K = jm = (mi, ...,nik) G : Ym=i m i = l} an< ^ d a is defined in (2.6). 



min {T a (g) : q € Plan [x, b]} = min ^d a mi ^ x n 



Proof. Since d a is a metric, it is easy to see that d a ^X)»=i m i$xn bj is a continuous function 

in (mi, • • • , mfe) on the compact subset K of Thus, the function d, 
achieves its minimum value in K at some rh* = (m*, ...,mt). That is, 



d a (a*,b) = min d a rm5 Xi ,h , 



(3.2) 



where a* = ^^m*^. Now, let G* be an optimal transport path in Path (a* ,h) with 
M a (G*) = d a (a*,h). By Xia |24| Lemma 7.1], G* has at least one compatible plan 
q* G Plan (a*, b) C Plan [x, b] and thus 

T Q (q*) < M a (G*) = d a (a*,b) . (3.3) 

For any q G Plan [x, b] , we have 



T a (g) > d a (a (g) ,b),by@ and §2A§ 

> rf„(a*,b),byj§ 

> T a (g*),by@. 
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This shows that 

T Q (q*) = min {T a (q) : q G Plan [x,b]} and T a (g*) = d a (a*,b) . (3.4) 



Thus, is a solution to the ramified optimal allocation problem (2.18). □ 

One implication of the above proposition is that there exists a close relationship be- 
tween optimal allocation plans and underlying transport paths. To further characterize the 
properties of an optimal allocation plan, we introduce the concept of allocation paths as 
follows: 

Definition 3.1. An allocation path from x to b is a transport path G G Path (a, h) for 
some atomic probability measure a supported on x. Denote Path [x, b] as the set of all 
allocation paths from x to b. An allocation path G* G Path [x, b] is optimal if 

M a (G*) < M a (G) for any G G Path [x, b] . 

By the definition of d a , equation ( |3.1[ ) can be alternatively written as 

min {T a (q) : q G Plan [x, b]} = min {M a (G) : G G Path [x, b]} , (3.5) 

which shows that the ramified optimal allocation problem corresponds to a problem of 
finding an optimal allocation path. The following lemma establishes the connection between 
optimal allocation plans and paths via the notion of compatibility. 

Lemma 3.1. Suppose G G Path[x, h], q G Plan[x,h] and (G,q) is compatible. Then, G 
is an optimal allocation path if and only if q is an optimal allocation plan with T a (q) = 
M a (G). 



Proof. If G is an optimal allocation path, then by (3.5), 

M a (G) = min {T a (q) : q G Plan [x, b] } . 

Since (G,q) is compatible, T a (q) < M a (G). Thus, 

T Q (q) = min {T a (q) : q G Plan [x, b]} = M a (G) , 

and q is an optimal allocation plan. 

Suppose q is an optimal allocation plan with T a (q) = M a (G), then 

M a (G) = T a (q) = min {T a (q) : q G Plan [x, b]} . 



By (3.5), G is an optimal allocation path. □ 

The next lemma presents a key property of an optimal allocation path. 

Lemma 3.2. Let G G Path [x, b] be an optimal allocation path from x to b. Then, for any 
i ^ s G {1, • • • , k}, Xi and x s do not belong to the same connected component of G. 

Proof. Assume and x s belong to the same connected component of an optimal allocation 
path G = {V(G),E(G),w : E(G) — > (0, +00)}, then there exists a polyhedra curve 7 
supported on G from Xi to x s . We may list edges of 7 as 

{eiei, • • • , e n e n } with e« = ±1 and et G E (G) . 

Here, Sj = 1 (or —1) if has the same (or opposite) direction as 7. Let 

A = min w (ei) > 0, 

Ki<n 
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and consider Gt := G + tj for t = ±A. Note that Gt is still in Path [x, b], and 

n 

M a (G t ) - M a (G) = [(«> (Ci) + tSiT - (w (ei)) a \ length (e t ) . 
i=l 

Thus, 

M Q (G A ) + M a (G_a) - 2M a (G) 

n 

= E t( w ( e *) + A ) a + ( w ( e *) - A ) a - 2 ( w ( e ^)) a ] Zen ^ ( e *) • 
i=i 

When a £ (0, 1), by the strict concavity of x a , we have 

( w ( e .) + A ) a + («, ( ei ) -X) a -2 (w ( ei )) a < 0. 

When a = 0, 

n 

[(«; (ei) + A) Q + (w (et) - X) a - 2 (w ( ei )) a ] length (e t ) 

= — {length (ei) : w (ei) = A} < 0. 
Thus, when a € [0, 1), we have 

M a (G A ) + M a (G_a) - 2M a (G) < 0. 

i.e. 

min {M a (Ga) , M a (G_ A )} < M a (G) , 
a contradiction with the optimality of G in Path [x, b] . □ 

The above lemma says that no two factories will be connected on any optimal allocation 
path. Alternatively speaking, on an optimal allocation path, each single household will 
receive her commodity from only one factory, i.e., each household is assigned to one factory. 
This result is attributed to the transport economy of scale underlying ramified transporta- 
tion technology. As seen in Section 2, an a £ [0, 1) implies the existence of transport 
economy of scale with transporting in groups being more cost efficient than transporting 
separately. Any allocation path on which some single household receives commodity from 
two factories can not be optimal because the planner would be able to reduce transportation 
cost by transferring production of one factory to the other. This transfer makes the benefit 
of transport economy of scale more likely to be realized as commodity for this household is 
transported in a larger scale on the path. 

The result that each household is assigned to one factory on an optimal allocation path 
motivates the following notion of assignment map. 

Definition 3.2. An assignment map is a function S : {1, ■ ■ ■ ,£} —> {1, ■ ■ ■ , k}. Let 

Map [£, k] be the set of all assignment maps. For any assignment map S G Map [£, k] 
and a 6 [0, 1), define 

k 

B a (5;x,b) :=J2d a (a^bj, 

i=i 
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where d a is the metric defined in (2.6) 



^ rij\S x .andbi= ^ njS y .. ( 3 -6) 

An assignment map S* G Map [£, k] is optimal if 

E Q (S*; x, b) < Eq, (S; x, b) for any S G Map [£, k] . 



Using the concept of assignment maps, Lemma |3.2| provides a partition result for an 
optimal allocation path. 

Proposition 3.2. Let G be an optimal allocation path from x to b. Then, there exists an 
assignment map S G Map [£, k] such that 

E Q (5;x,b) = M a (G). 

Moreover, for a^ and hi given in {3.6^ , the path G can be decomposed into the sum of k 
pairwise disjoint transport paths 

k k 
G = ^Gi with M Q (G) = M a (Gi) , (3.7) 

i=l i=l 

where each Gi G Pai/i (aj, bj) is an optimal transport path with M a (Gi) = d a (a^, bj . 

Proof. Let G be an optimal allocation path from x to b. By Lemma [3.2[ each Xi determines 
a connected component Gi of G. Thus, 

k k 

G = Y,Gi with M Q (G) = Y M * (Gi) • (3.8) 

i=l i=l 



For each j G {1, • • • ,£}, yj is clearly connected to some Xi on G. By Lemma 3.2, such an 
Xi must be unique for each j. Thus, we may define 

S (j) = i, if Uj is connected to x$ on G. 

This defines an assignment map S G Map [£, ft]. Also note that Gi is a transport path from 



aj to hi, where a^, hi are given in (3.6). Since G is optimal as in (3.5), each Gi G Path (a^, hi) 
must also be optimal, and hence 

M a (Gi) = d a (a;,bj). 

Therefore, 

k k 

M a (G) = M « ( G = E ^ ( ai ' bi ) = E - ^ x ' b ) ' 

i=i i=i 

□ 

Figure [4] illustrates a partition result under an optimal assignment map. This spatial 
economy consists of factories {x±, x%} and households {y±, • • • , 2/145} with equal demand nj = 
The figure shows a clear partition in both the allocation path and households: First, 
the allocation path is decomposed into two disjoint sub-transport paths originating from 
factory x\ and X2 respectively. Second, the 145 households are divided into two unconnected 
populations centered around the factory from which they receive the commodity. 
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5.5 6 6.5 



7.5 e 3.5 



Figure 4. A partition in households b and allocation path G. 
Note that any assignment map 5 £ Map [£, k] provides a partition among households 



b =E b * 



i=i 



where all households in bj given in (3.6) are assigned to a single factory located at Xj. For 



each i, let G{ £ Path (aj, bj be an optimal transport path which transports the commodity 
produced at factory i to households in bj. This yields an allocation path from x to b 

k 

G s = J2 Gie Path I x ' b l ' ( 3 - 9 ) 



i=l 



which is compatible with the allocation plan q$ given in (2.15). Thus, 

k k 



B a (S; x, b) = < (ai, bj = ^M U > M a (G s ) > T a (q s ) . (3.10) 



i=l 



i=l 



By Proposition 3.2, any optimal allocation path is in the form of (3.9) with respect to its 



associated assignment map. 

Now, we state the main results of this section as follows: 

Theorem 1. Given a subset x = {x\, X2,~. Xk} in X, an atomic probability measure h as 
in (2.1), and a parameter a 6 [0, 1). 

(1) An allocation plan q £ Plan [x, b] is optimal if and only if there exists an optimal 
assignment map S £ Map [£, k] such that q = qs- 

(2) An allocation path G £ Path [x, b] is optimal if and only if there exists an optimal 



assignment map S £ Map [I, k] such that G = Gs for some Gs of S defined in (3.9). 
(3) Moreover, 

min T a (q)= min E a (S';x,b)= min M a (G). (3.11) 

q£Plan[xL,b] SeMap[£,k] G&Path[yL,b] 
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Proof. We first show the equivalence in ([!]): 

"^=^" . Suppose q G Plan [x, b] is an optimal allocation plan. Let G G Path [x, b] be 
an allocation path that is compatible with q and M a (G) = T Q (q) . Since q is an optimal 
allocation plan, by Lemma |3.1[ G is an optimal allocation path. By Proposition |3.2[ there 
exists an assignment map S such that 

E Q (S; x, b) = M a (G) = T a (q) . (3.12) 

Since (G, q) is compatible, by ( 2.12| ), we may express 

k t 

G = J2J2 qi i " gi r 



i=i i=i 



For each j, by Lemma 



3.2 



we have g s j = whenever s ^ S(j). By the compatibility of 
(G, q), we have g S j = whenever s ^ S (j). Thus, for i = S (j) , 



^2 Qsj = nj, for each j. 



This shows that q = qs- 

For any S G Map [£, fc], by j3.10| ) and ( |3.12| , 

E Q (5; x, b) > T a (q § ) > T Q (q) = E Q (S; x, b) . 

Thus, S is an optimal assignment map. Since E Q (S';x,b) = T a (q), we have the first 
equality in (3.11 ). 

"<£^". Suppose S G Map [£, k] is an optimal assignment map. Then, by (3.10) and the 
first equality in ( 3.11| ), 

T a (q s ) < E Q (5; x, b) = min {E a (§; x, b) : S G Map [i, fc]} 

= min {T a (g) : g G PZan [x, b]} . 

Therefore, 05 G PZan [x, b] is an optimal allocation plan. This completes the proof of ([!]). 
The second equality of (3.11) follows from the first equality of (3.11) and (3.5). 
We now show the equivalence in (|2j): 

"^=" Suppose G = Gs for some optimal assignment map S G Map [£, k] , then 

k k 

Mq, (Gg) < ^M a (G J ) = ^d a (a i ,b J ) = E a (5;x,b) 



i=l 



min E a ( S; x, b 



min M a (G) 

GePat/i[x,b] 



by the optimality of S and (3.11). Thus, G = G5 is an optimal allocation path. 



"^=^" Suppose G G Pai/t [x, b] is an optimal allocation path. By Proposition 3.2 G = Gs 
for some S G Map [£, k] with E a (S; x, b) = M a (G). Then, the optimality of S follows from 
the optimality of G and ( |3.11 ). □ 

Theorem [T] shows that in the ramified optimal allocation problem, there exists an equiv- 
alence between optimal allocation plan and optimal assignment map. This result has an 
analogous counterpart in Monge-Kantorovich problems, but has not been observed in cur- 
rent literature on ramified transport problems. An implication of this theorem is that one 
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can instead search for an optimal assignment map in order to find an optimal allocation 
plan. Each optimal assignment map S G Map [£, k] would give an optimal allocation plan 
qs G Plan [x, b] as in (2.15). For this consideration, the rest of this paper will focus atten- 
tion on characterizing the various properties of optimal assignment maps. 

4. Properties of Optimal Assignment Maps via Marginal Analysis 

In this section, we develop a method of marginal transportation analysis and use it to 
study the properties of optimal assignment maps. 

We first formalize a concept of marginal transportation cost for a single-source transport 
system. Let G G Path (m (c) 5o, c) be a transport path from a single source O to an atomic 
measure c of mass m (c). For any point p on the support of G, we set 

w(e), if p is in the interior of some edge e G E (G) 

0(p):=im(c), ifp = , (4.1) 

,E eeE(G ), e+=P ™(e), if P eV(G)\{0} 

which represents the mass flowing through p, where e + denotes the ending endpoint of edge 
e G E (G). Since G has a single source and contains no cycles, for any point p on G, there 
exists a unique polyhedral curve 7 P on G from O to p. Moreover, for any point s G 7 P , it 
holds that 

m(c) >0(s) >0(p). (4.2) 

When the mass at p changes by an amount Am with Am > — 9 (p) , then the mass flowing 
through each point of 7 P also changes by Am. As a result, the corresponding increment of 
the transportation cost is 

AC G (p,Am): = M a (G + (Am) j p ) - M a (G) (4.3) 
= f {9(s) + Am) a -(9(s)) a ds. 

The marginal transportation cost at p via G is defined by 

MC G (p):= lim AC f Am) = J (9(s)Y-Us. 

The following proposition establishes some properties of the function ACq (p, Am). Those 
properties are the key elements of marginal transportation analysis used to study optimal 
assignment maps later. 

Proposition 4.1. For any G G Path (m (c) 5o, c) and p on G, we have 

AC G (p, -Am) = -ACq (p, Am) , for Am G [-9 (p) , 9 (p)] , (4.4) 
where G = G — (Am) 7 P . Moreover, for any Am > 0, we have 

[(m (c) + Am)" - m (c) a ] J ds < AC G (p, Am) < [(0 (p) + Am)° - 9 (p) a ] ! ds. (4.5) 

•>1p Jip 

If, in addition, G G Path (m (c) 5o, c) is optimal in ^2Jm , then 

d a ((m (c) + Am) 5 G , c+ (Am) <5 p ) - d Q (m (c) 5 , c) < AC G (p, Am) (4.6) 
/or any Am > — 9 (p) with p on G. 
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Proof. Let G = G — (Am) j p , then when Am G [—9 (p) , 9 (p)], 
AC G {p,-Am) = I (9 (s) - Am) a - (9 (s)) a ds 

= - [ [((9(s)-Am) + Am) a -(9(s)-Am) a ]ds 
= -AC G (p,Am). 

For any Am > 0, note that the function / (t) := (t + Am)° — t a is monotonic non-increasing 
on t > when < a < 1. By (4.2), we have 

(m (c) + Am)" - (m (c)) a < (9 (s) + Am) Q - (0 (s)) Q < {9 (p) + Am) a - (0 (p)) Q . 

Thus, by p~3| , inequalities (j45j) hold. 

When G is also optimal, we have M Q (G) = d a (m (c) So, c). For any Am > — 9 (p), since 
G + (Am) 7 P G Pat/i ((m (c) + Am) So, c+ (Am) S p ) , we have 

AC G (p,Am) = M a (G + (Am) 7p )-M Q (G) 

> d a ((m (c) + Am) So, c+ (Am) 5 P ) — (i a (m (c) So, c) . 

□ 

Now, we apply marginal transportation analysis developed above to study properties 
of an optimal assignment map. Let G G Path [x, b] be any optimal allocation path. By 
Theorem [TJ G must be in the form of (3.9), which is simply a disjoint union of single-source 
paths GVs. For any p on G, there exists a unique i such that p is on G{. Thus, we can 
define the corresponding 9i (p) and ACq 1 (p, Am) as in (4.1) and (4.3). Then, we set 

9 (p) := 9i (p) and AC G (p, Am) := AC Gi (p, Am) 

for any p on the support of G. 

For any a G (0, 1) and a > e > 0, define 

Pa (o-,e):=(°Y • (4-7) 

The function p Q (a, e) is decreasing in a and increasing in e; also < p a (a, e) < 1 for any 
o > e > 0. For a = 0, set (c, e) = when a > e > and po ("" 5 e ) = 1 when cr = e > 0. 



Proposition 4.2. Suppose G G Pai/i [x, b] is an optimal allocation path as given in (3.1). 
Let p be a point on G s for some s G {1, • • • , k} and Am G (0, 9 (p)]. For any p* on G with 
7p H 7 P * having zero length, we have 

AC G {p,-Am) + (Am)" \\p — p*\\ + AC G (p* , Am) > (4.8) 

and 

\\P-P*\\ +Pa (9(p*) + Am, Am) / ds > p a (m (b s ) , Am) / ds, (4.9) 
where \\-\\ stands for the standard norm on W m . In particular, for any i G {1, • • • , k}, 

\\p - Xi\\ > p a (m (b s ) , Am) ds > p a (m (b s ) , Am) \\p - x s \\ . (4.10) 



Moreover, suppose p* is on G{ with i ^ s and Am < 9 (p*), then 
n * u , Pa(9{p*) + Am,Am) 

\\p-p H / /u \ a — \ — \\P ~ Xi\\> p Q (m(b s ) , Am) \\p- x s \\ . (4.11) 

p a (m (hi) , Am) 
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Proof. Let G = G— (Am) 7 p +(Am) [p,p*] + (Am) 7 P », where [p,p*] denotes the line segment 
from p to p*. Then, when the intersection of polyhedral curves 7 P n j p * has length zero, we 
have 



AC G (p,-Am) + (Am) Q ||p-p*|| + AC G (p*,Am) 
f [{9 (a) - Am)" - (s) Q ] (is + (Am)" ||p - p*|| 



+ / |(0(s) + Am) a -0(s) a ]ds 
> M a (g) - M a (G) > 0, by the optimality of G. 



To prove (|4.9|), we observe that 

\\p — p* \\ + p a (9 (p*) + Am, Am) / 



(Am) c 



(Am) a ||p - p*|| + [(9 (p*) + Am) a - (9 (p*)) a ] I 



ds 



> (Am)° l|p-p*[[ + ACq (p*,Am) , by (|4~5 

> —ACq (p, -Am) , by g 



• 5) 



= AC 6 (p,Am), by Q 
> [(m (b s )) a - (m (b s ) - Am)"] / ds, by gj 

= (Am) a p Q (m(b s ),Am) / ds > (Am) a p a (m (b s ) , Am) \\p — x s 

J 7n 



'7p 

Xi for some i, (4.9) becomes (4.10) as f^^ds = 0. 



In particular, when p* 
Now, suppose p* is on Gi with i ^ s. Apply (4.10) to p*, we have 

Hp* - xj 



p a (m(bi) , Am) 



> 



ds. 



Therefore, 



p a (0(p*)+Am,Am) , 

p — p H 7 — 77—, — 7 — \ — Hp - Xi 



p a (m (bi) , Am) 

> ||p - p* || + p a (9 (p*) + Am, Am) / 

> p a (m (b s ) , Am) / ds, by 

Ap 

> p a (m(b s ) , Am) ||p - x s \\ ■ 



ds 



□ 



Proposition 4.2 rests on the key inequality (4.8), which follows intuitively by standard 
marginal argument. For convenience of illustration, let p denote the location yj of household 
j who is connected to factory s by some curve j p . A planner will find it not optimal to choose 
an allocation path G such that AGg (yj, —Am) + (Am)™ \\yj — p*\\ + AGg (p*, Am) < 
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for p* on some Gj. It is because in this case the planner has a less costly alternative by 
transferring Am amount of production from factory s to factory i and transporting this 
additional Am units of commodity from factory i first to a stopover point p* via curve "y p * 
and then directly from p* to household j. It's clear that this strategy will send the same 
amount of commodity to household j as before. However, by the inequality, the reduc- 
tion in transportation cost —ACq (p, —Am) on curve 7 P exceeds its increase counterpart 
ACg (p* , Am) + (Am)" \\yj — p*\\ , which cannot be the case for an optimal allocation path. 

By means of this proposition, we easily obtain the following results regarding the prop- 
erties of optimal assignment maps. 

Theorem 2. Suppose S G Map [£, k] is an optimal assignment map. Let j G {1, • • • , £} and 
s € {!,-•• ,k}. If 



Vj G F's(nj) := Iz G 
then S (j) ^ s. If 



: mm \\z — x. 



ill < p a (m(b s ) ,Uj) \\z - X s || j , 



Vj G ( n j) 



z G 



x s \\ < min p a (m(bj) ,nj) 



(4.12) 



(4.13) 



then S (j) 



s. 



Proof. Assume yj G Fg (nj) but S (j) = s. Let G = Gs be an allocation path as in (3.9). 
By Theorem [TJ G G P ath [x, b] is an optimal allocation path. Clearly, yj is on G s when 



S (j) = s . Apply (|4. 10 ) to p = yj and Am 

mm \\yj - Xi\\ > pa 
ie{l,...,fc} 



nj, we have 



> (b s 



which contradicts (4.12). 



On the other hand, for each i ^ s, if yj G Q% i n j)i then 



mm \ \yj 



< 



Xq 



< Pa (TO (hi 



Wj 



By (4.12), yj G F 



and S (j) 7^ i for such i ^ s. Thus, S (j) = s. 



□ 



Intuitively speaking, inequality ( |4.12 ) says that if household j locates "closer" to some 
factory than factory s, then the planner will not assign her to the factory s. Here the 
relative closeness is weighted by a number p a (m (b s ) , nj). When the production m (b s ) at 
factory s is low, due to the transport economy of scale, one would expect that the planner 
would less likely assign household j to factory s. This predication is justified by Theorem 
[2] because in this case, inequality (4.12) becomes more likely to hold as p a (m(b s ) , nj) is 
high. The later part (4.13 ) of the theorem states a special case that if household j is located 



uniformly closer to a factory s than to other factories, then she will be assigned to factory 
s under any optimal assignment map. 

We now give a geometric description of the sets F s s (nj) and Q s s (nj). 



Lemma 4.1. For any constant C G (0, 1), the set 



{x G M m : \\x - Xi\\ <C\\x- x s \\} = B | r, - 
where B (x, r) denotes the open ball {z G 



C 2 



(J2 ( X * ' 
II < r }- 



c 



Xq 



1-C 2 
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Proof. Indeed, 



| x x% 1 1 C 1 1 x x g 1 1 

i H 2 . n2 ii I ■ 2 

•i - ^ " y Ju Ju g 

x 1 1 2x * ~\~ 1 1 x^ 1 1 -^C C [ 1 1 x 1 1 2x * x g ~\~ 1 1 x 



2x 



1-C 2 

c 2 



< 



C llx.. 



1-C 2 

2 ^ 2 



Clearly, the ball B ( Xi + 



1-C 2 



< 



C 2 



;i - c 2 y 



Xo X a 



□ 



IX?; Xc 



c 



.» ,. — >\\ | given above contains xi but not 

x s . By Lemma 4.1 (4.12) says that geometrically, if household j lies in the union of k — 1 
open balls 

\2 



W 



\Xi X s ) , 



w 



1 - (w s j) 



\Xa Xj\ 



where w s j = p a (m(b s ) ,nj), then S (j) ^ s. Also, (4.13) says that if the household j lies 
in the intersection of (k — 1) balls 



n s s (n 3 ) = f]Blx s + - 



11/ 



{wij)' 



(x s Xi) 



W; 



i - (wtjY 



Xi Xg 



(4.14) 



for some s, then S (j) = s. 

Since m(b s ) < m (b) = 1 and the function p a (-,nj) is decreasing, we have p a (l,n.j) < 
p Q (m (h s ) , rij) and thus 



F s (rij) : = < z : min\\z — Xi\\ < p a (l,rij) \\z — x s \\ > C. Tg (rij) , (4-15) 

I i+s J 

fl s (rij) : = < z : \\z — x s \\ < p a (1, rij) min \\z — Xi\\ > C U s s (nj) (4-16) 

for any optimal assignment map S. Note that the set F s (rij) (or Q s (rij)) is still the union 
(or intersection) of k — 1 balls in M m , and is independent of S. For example, the sets 
Sl s (rij) with rij = 0.8 and re.,- = 0.5 are given by Figure [5j where x\ = (0, 0), X2 = (2, 0) and 
x 3 = (l,2). 

Corollary 4.1. For any optimal assignment map S G Map [i, k] and s G {1, • • • ,k}, if 
Dj G F s (rij), then S (j) ^ s. If yj G Q s (nj), then S (j) = s. 



Proof. It follows from Theorem [2j ( 4.15| ) and (4.16). 



□ 



This corollary shows that if the household j falls into the region £l s (nj) of some factory 
s, then she will be assigned to this factory under any optimal assignment map S. As a 
result, if all households belong to the union of regions f2 s (nj) of factories s except factory 
i, then factory i will not be used. Note that as p a (l,re?) is increasing in nj, the size of the 
region fl s (nj) increases with nj as shown in Figure [5j 



20 



QINGLAN XIA AND SHAOFENG XU 




Figure 5. An example of the sets Q s (rij) with rij = 0.8 (blue) and nj = 0.5 
(red) when a = 1/2. 



Theorem 3. Suppose S G Map[£,k] is an optimal assignment map, h and j G {1, 
with rij <rih- If S (h) = s* ^ s for some s G {1, • • • ,k} and 



z G 



Phil + 



< p a (m (b a ) , nj 



p a {n h +nj,nj) 



\\yh 



X s * 



Z - X, 



(4.17) 



then S (j) ^ s. If S (h) = s for some sG {!,■•■ , k} and 



Vj : Of (rij) := { 



z G 



\\ z - vh\\ + p p}(S n J',n-) \\y h -s\ 

< p a (m(bi) ,n,j) \\z - Xi\\ 



(4.18) 



then S (j) = s. 



Proof In the first scenario, assume S (h) = s* 7^ s, yj G F S g h (rij) but S (j) = s. Then 



yh\\ + — / / U \ — a ll^fc ~ x * 



> 



Vh\\ + 



Pa (m (b s ») ,nj) 

p a (g(y ? ) + n 3 -,n 3 -) as n „ 

p Q (m(b s *) ,rij) 



> p a (m(b s ) ,n,) ||yj - x s || , by (4.11), 



a contradiction with y 3 G Fg (n-,-). Thus, S (j) 7^ s. 
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Now, in the second scenario, assume S (h) = s, yj G Qg h (rij) but S (j) = i* for some 
i* 7^ s. Then, 

WVj -yh\\-\ t — 77-\ r hh ~ x s \\ 

p a (m(b s ) ,Uj) 

> hj - VhW + Pa ^7^7 ' + nj '" J ^ || y/l - a; s || as > n h , 

p a (m (b a ) , rij) 

> p a (m (bj») , rij) hi - Xj* || , by fl4.1ip , 

a contradiction with G Og (rij) as m (bj*) > rij and i* 7^ s. Thus, S (j) = s. □ 

The first part of Theorem [3] says that if some household h is not assigned to factory s, 
then any other nearby household j (i.e. within the neighborhood region Fg h (rij)) with 
a smaller demand will also not be assigned to factory s. The second part says that if 
some household h is assigned to factory s, then any other nearby household j (i.e. within 
the neighborhood region Q s ^ h (rij)) with a smaller demand will also be assigned to factory 
s. These findings agree with the intuition that grouping with nearby households of large 
demand would make it more likely to realize the benefit of transport economy of scale. 

5. Properties of Optimal Assignment Maps via Projectional Analysis 

As seen in the previous section, under an optimal assignment map, a household will be 
assigned to some factory if she lives close to the factory (Theorem [2]) or she has some nearby 
neighbors assigned to the factory (Theorem [3]). In this section, we will show a reverse result 
(Theorem [4| using a method of projectional analysis. 

Throughout this section, we consider the projection map from M m to the line {p + tv : t € M 
for fixed points p,v £ M. m with \\v\\ = 1. Under this map, each point z 6 IR m is mapped to 
p + 7r (z) v with 

tt(z) = (z-p,v) , (5.1) 
where (•, •) stands for the standard inner product in W 1 . For instance, when p = (0, • • ■ , 0) 6 
M. m , and v = (1,0,- •• ,0) € M m , for each z = (z\,--- ,z m ), tt (z) = Z\ gives the first 
coordinate of z. 

We start with two lemmas regarding properties of a single-source transport system. These 
lemmas will play a crucial role in establishing Theorem [4] later. 

Lemma 5.1. Suppose 

c = c j^Zj, Cj > for j in a finite set (5-2) 
jee 

is an atomic measure on X C W 71 , and P,Q £ X . If there exists a t\ such that 

7r (P) > t\ > max7r (zj) or tt (P) < t\ < min7r (zj) , 
jee jee 

then 

d a (c, m (c) 5 P ) - d a (c, m (c) 5 Q ) > m (c) a (|tt (P) - h\ - CR - |tt (Q) - h\) , 

where 

C = 2 i_ (m _i)(i_ Q ) _ 1 + !> ( 5 - 3 ) 

and 

R = m&x{\\z-p-Tr(z)v\\ : z G {zj : j G Q}U{P,Q}} . (5.4) 
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Proof. Without loss of generality, we assume that ir (P) > t\ > max^ge 7r (zj). Let G G 
Pat/t (c, itl (c) dp) be an optimal transport path. As in (2.11), there exists a unique curve 
7 2j . from P to Zj for each j. Since it (P) > ti > it (zj), there exists a point vj on the curve 
7 2j such that 7r (wj) = t\. Let 

C = Cjfi^ , 

iee 

then, by the optimality of G, we have 

d a (c, m (c) 5p) = d a (c, c) + d a (c, m (c) 6p) . (5.5) 

For 0\ = p + tiw and P given in ( |5.4[ ), we consider the set 

B R (Ox) := {z G M m : tt (z) = ti, ||z - Oi|| < P} , 

which is an (m — 1) dimensional ball perpendicular to the line passing through p in the 
direction v. By means of (5.4), the measure c is supported on Bp (Oi). Let 

P = P+(t 1 -TT(P))v 

be the projection point of P onto the ball Bp (0±). One can show as in Xia and Vershynina 
|31[ Lemma 2.3.1] that 

d Q (c,m (c) 5p) > d a (m (c) 5p,m (c) Sp) = m (c) a |vr (P) — ti|. 

On the other hand, by Xia [24[ Theorem 3.1], we have 

d a (c,m(c)5 Q ) < d a (c,m{c)6 0l ) + d a (m(c) 5 0l ,m{c) 5 Q ) 
< m(c) a CP + m(c) Q (P + K(Q)-t 1 |), 

where 

~ = y/m - 1 

- 2 l-(m-l)(l-«) _ 1' 

Therefore, 

d a (c, m (c) 8 P ) - d a (c, m (c) 8q) 

> d a (c, m (c) Sp) — [d a (c, c) + d a (c, m (c) Sq)] , by the triangle inequality 
= d a (c, m (c) 5 P ) - d a (c, m (c) Sq) , by (pT5|) 

> m (cf |tt (P) - ti| - m (c) a (CP + \ir (Q) - h\) . 

□ 



Lemma 5.2. Zei c be an atomic measure as given in (5.2) and O G X . If the set O is 

decomposed as the disjoint union of two nonempty subsets 

e = ©iiie 2 , 

then, for any optimal transport path G G Path (m (c) So, c), there exist a vertex point P G 
V (G) and a decomposition of each Oj/ 

Oi = Oj II <8>i with &i nonempty, i = 1,2 

such that G can be decomposed as 

G = G 1 + G 2 + G Z (5.6) 
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where for i = 1,2, Gi is an optimal transport path from m (£j) 5p to Cj for 

G3 is an optimal transport path from m (c) 5o to c+ (m (£1) + m (£2)) Sp for 

jeOiue 2 

and {Gi}f =1 are pairwise disjoint except at P. Moreover, (5.6) implies 

M a (G) = M a (G 1 ) + M a (G 2 ) + M a (G 3 ) (5.8) 

and by the optimality of G, it follows 

d a (m (c) <5o, c) = d a (m (£1) <5p,£i) + d a (m(£ 2 ) 5p,c 2 ) (5.9) 

+d a (m (c) 5 , c+ (m (£1) + m (£ 2 )) 5 P ) . 

Proof. For any z on the support of G, since G is a transport path from a single source O, 
there exists a unique curve "y z on G from O to z. Also, it is easily observed that 

if z lies on j z for some z, then j z is the part of 7 2 from O to z. (5.10) 

Now, let Tj be the union of all curves 7^. with j £ 0, for i = 1,2, and set 

r = r!nr 2 . 



By (5.10), if z G V, then j z C T. This shows that T is a connected subset of the support of 
G containing O. Since T contains no cycles, it is a contractible set containing O. Then, by 
calculating the Euler characteristic number of T, we have either T = {O} or T has at least 
two endpoints (i.e. vertices of degree 1). 

If T = {O}, then set P = O, and 6; = 6, for i = 1,2. If T ^ {O}, pick P to be an 
endpoint of T with P ^ O. Since P G T C Tj, the set 

6i := { j S Oj : P G 7,.} / 0, for i = 1,2. 

For any j G G>i, P divides the curve 7^. into two parts: 7^ from O to P and 7^ from 
P to Zj. Since P is an endpoint of T, we have 

; 7 g\{P})nr = 0. 



For i = 1, 2, define £^ using (5.7) and denote the part of G from m (£j) <5p to £j by G^. The 
rest of G is denoted by G3 = G — (Gi + G 2 ). Then, by construction, {G{\\ =1 are pairwise 
disjoint except at P, and thus (5.8) holds. By the optimality of G, each Gi must also be 
optimal for i = 1,2, 3, which yields (5.9). □ 

The following theorem states that: under an optimal assignment map, a household will 
be assigned to some factory only when either she lives close to the factory or she has some 
nearby neighbors assigned to the factory. In the first situation, the planner takes advantage 
of relative spatial locations between households and factories; while in the second situation, 
the planner takes advantage of group transportation due to transport economy of scale 
embedded in ramified transport technology. 

Let S G Map [£, k] be an optimal assignment map. For each i G {1, • • • , k}, define 

Vi = {y j :S{j)=i}U{x i } (5.11) 
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and 

Ri = max{||z - p - ir (z) v\\ : z G {yj : S (j) = i} U {x\, ■ ■ ■ ,x k }} ■ (5.12) 

Theorem 4. Suppose S G Map [£, k] is an optimal assignment map. Then, for any i G 
{1, • • • , k} and j G {1, • • • ,£} with S (j) = i, there exists z G {yj}, such that 

< \w ( yj ) - ir (z)\ < 2CRi + mm\ir(xi*) - tt ( Vj )\ (5.13) 



and 7r (z) is between tt (yj) and tt (xi), where C and Ri are the constants given in (5.3) and 
(5.12) respectively. 

Proof. Without loss of generality, we may assume that tt (yj) < tt (xi). Let 

|vr (xi*) - it (yj)\ = min |vr (xj»») - vr (yj)\ 



for some i* ^ i. We want to prove (5.13) by contradiction. Assume for any z G ^i\{yj} 
with TT (yj) < tt (z) < TT (Xi), 

\tt ( yj ) - TT (z)\ > 2CRi + \<IT (Xi*) - TT (yj)\ , 

i.e. 

vr (z) - tt ( Vj ) > 2CR l + \tt (x?) - tt (yj)\ . 
Then, there exists a real number t<± such that 

TT (z) > t 2 > tt (yj) + 2CR { + \tt ( Xi *) - TT ( V j)\ (5.14) 
whenever tt (yj) < tt (z) < tt (xi). In particular, 

rr(xi) > t 2 > tt ( yj ) + 2CRi + \tt ( Xi *) - tt ( Vj )\ (5.15) 
as Xi G As a result, S^ 1 (i) can be expressed as the disjoint union of two sets: 

sr 1 ® = eine 2 , 

where 

6! := {h G S- 1 (i) : vr (y h ) < tt ( Vj )} , 9 2 := {h G S' 1 (i) : vr (y h ) > t 2 } . (5.16) 
Clearly, j G @i. If 9 2 = 0, then S' 1 (i) = 6i and thus 

vr (xi) > tt (yj) > max vr (y h ) . 

hG5 _1 (j) 



Let aj and bj be given as in (3.6). By Lemma 5.1 with t\ = tt (yj), we have 

d a (bi,m(bi)S Xi ) -d a (hi,m(bi)5 Xit ) 

> m(bi) a (\TT(x l )-TT(y J )\-CR i -\TT(x i *)-TT(yj)\) 

> m (bi) a (tt (x^ - tt (yj) - CRi - |vr (x^) - tt (yj) \) 

> m (bi) a CRi > 0, by ( foTl5| ), 

a contradiction to the optimality of S. Thus, G 2 ^ 0. 

Let Gi G Path (a^, bj) be an optimal transport path. Then, by Theorem [T] and optimality 
of S, G = Y2i Gi is an optimal allocation path. Since both 0i and 2 are nonempty, by 



setting B = 5 1 (i) = Gi II @ 2 ,0 = Xi and c = bj in Lemma 5.2, there exists a point 
P G V (Gi) such that Gj can be decomposed as 

Gi = G^ + GP + GP 
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with 



M a (Gi) = M a (g, (1) ) + M a (of >) + M Q (g\ 



(3) 



(5.17) 



Here, for h = 1,2, G-^ is an optimal transport path from m [h\"'' ) <5p to b^"'' for some 
positive atomic measures b^ with spt ^b^^ C : /i £ fc | and 



:(*) 



:C0 



Gf } e Path I m (bi) ff^, hi - I b^ + hY> ) + ( m ( b^ ) + m ( b^ ) ) 5 P 



(l)^K( 2 ) 



;(2) 



If 7T (P) > t2 — CRi, then we can modify G into another allocation path G by just 
replacing the corresponding transport path from factory i to households b^ with an optimal 
transport path from factory i* to b\ ■ More precisely, we replace Gj by 



(5.18) 



where G\ is an optimal transport path from m ( b-^ ) 5 X t to b f w , Gf ! = Gf' and 

Gf =Gf>- m (m 1P , 



-(1) A(2) _ n (2) 



(5.19) 



where 7p is the curve on G from Xi to P. Equation (5.18) and (5.19) imply respectively 



M a (Gi) < M Q ( G- 1 -' ) + M Q ( Gf ) + M a ( G 



,(2) 



*(3) 



and 



M„ ( G? } ) < M„ ( G 



,(3) 



Consequently, by ( 5.17[ ), 

M Q (d) -M a (G 



> M a [ GP )-M a ( GY> ) + ( M a I G)°> )-M a (G 



,(3) 



V(3) 



> M„ G 



M, 



da [ml bf ] ) dp, by ) - d a I m ( b^ ) <k„ , b 



(5.20) 



where the last equality follows from the optimality of both G^ and G^ . 
Since a < 1, for G = ^» + we nave 



M Q (GJ < M « (G.) + M a (Gi) . 



Due to the optimality of G, equation (|3.8l) says 



M Q (G) = ^M Q (G s ) + M a (Gi 
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As a result, 

M a (G) - M a (g) 



> M Q (Gi) — M a I G 



> d a (m (bf } ) 5 P , b 4 (1) ) - d a (m (hf ] ) 6 Xi * , b«j , by 

> m (b? } )° (|vr (P) - vr (%■) | - Ci^ - |vr (x;*) - vr (%•) |) , by Lemma [ju] 

> m (b^Y (vr (P) - vr (jfe) - CRi - | vr (a*) - vr (%■) |) 

> m (b[ l) y (t 2 - it fa) - 2CRi - \ir (xr) - it ( Vj ) |) > 0, by ( pi] ). 

Thus, M a (G) > M a ^G^ , which contradicts the optimality of G, and thus the inequality 
(5.13) must hold. 

If it (P) < t 2 — CRi, then let Q be the first point of jp with it (Q) = t 2 . We can modify G 
into another allocation path G by just replacing the corresponding transport path from the 

~ (2) ~ (2) 

point Q to households h\ with an optimal transport path from Q to b- . More precisely, 
we replace Gi by 

where Gf } = G^, <5f } is an optimal transport path from m (b- 2 ^ 5q to b^ 2 \ and 

Gf) = G?)-m(bf) 7QPj 

where jqp is the part of the curve 7p from Q to P. Similar arguments as in the previous 
case show that 

M a (G) - M a (G) 

> d4mfbf))^,bf)-d4mfbf)j Q ,bf 

> m 



bf } J (i 2 - vr (P) - CPi) > 0, by Lemma |5.l|. 



Thus M Q (G) > M a (G) , which contradicts the optimality of G. Therefore, the inequality 
d5.13[ ) must hold. □ 

The following corollary states a scenario when a factory is located far away from the 
community of households, a planner will never assign any production to this factory under 
any optimal assignment map. 

Corollary 5.1. Suppose for some i £ {1, • • • , k}, 

\ir (xi) — it (yj) | > 2GP + min \tt (x^) — it (yj) | (5.21) 



for each j = 1, • • • , I, where C is the constant given in ]5.3[ ! and 

R = max{||z - p - it (z) v \\ : z € {y\, ■■■ ,ye}U {x 1 , ■ ■ ■ ,x k }} . (5.22) 
Then, S 1-1 (i) = for any optimal assignment map S £ Map [I, k\. 
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Proof. Assume there exists j G {1, • • • ,£} with S (j) = i. Without loss of generality, we 
may assume tt (xi) > tt (yj). Thus, there exists an j* G S -1 (i) such that 

tt (yj*) = max{7r (y h ) : tt (xi) > tt (y h ) , S (h) = i} . 

For this j* , by Theorem |4j there exists a z G with tt (yj*) < tt (z) < tt (xi) satisfying 
(5.13). By the maximality of n(yj*), z ^ yh for any y^ G Thus, z = xi and by (5.21), 

tt (xi) > tt (yj*) + 2CR + min \tt (xi*) — tt (yj*) \ 



> tt (yj* ) + 2CRi + min \tt (x^ ) - tt (yj* 

i*^i 



a contradiction with (5.13). 



□ 



The next corollary shows an "autarky" situation: if households and factories are located 
on two disjoint areas lying distant away from each other, then the demand of households 
will solely be satisfied from factories within the same area. 



Corollary 5.2. Suppose 

{tt (xi) , • • • , tt (x k ) , tt (yi) , • • • , tt (ye)} n (h,t 2 ) 
for some ti, t2 G M with 

t 2 = h + 2CR + a,a>0 



(5.23) 



where the constants C and R are given in (5.3) and (5.22). If 

{tt (x x ) ,tt (x k )} n (h - a, h] / and {tt (xi) , ■ ■ ■ , tt (x k )} n [t 2 ,h + <r) ^ 0, (5.24) 

then for any optimal assignment map S G Map[£, k] and interval I = (— oo,ti] or [£2,00), 
we have 

tt ( Vj ) G I ^> tt (x S (j)) G I, 

for j = !,-■■ ,£. 

Proof. It is sufficient to prove that if tt (xi) G [t 2 ,oo) for some i, then the set 

{tt (y h ) ■ S (h) = i, tt (y h ) < tx) 
must be empty. Indeed, if not, pick 

ir(yj) = max{Tr(y h ) : S (h) = i, tt (y h ) < tx} 
for some j. By Theorem |4j there exists z G ^/j such that tt (yj) < tt (z) < tt (xi), and 
tt (z) — tt (yj) < 2CRi + min \tt (xi*) - tt (yj)\ 

i*^i 

< 2CRj + tx-TT ( Vj ) + a, by gjg 

< t 2 -TT (yj) . 



Thus, tt (z) < t 2 . On the other hand, the maximality of tt (yj) and (5.23) yield tt (z) > t 2 , 
a contradiction. □ 



As a direct application of Corollary 5.2, the next corollary states that households living 
in a relatively isolated area are more likely to receive their commodity from local factories. 
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Corollary 5.3. Let t l < t 2 < < t% be real numbers with 

t~ = t~ + 2CR + a and tf = tf + 2CR + a, a > 0, 
where the constants C and R are given in fOU and (SJ§. If 

{vr(xi),-- - ,vr(x fc ),vr(yi),--- , vr (y e )} n ((tj , tj) U (#,#)) =0, 

{vr (xi) , • • • , vr (x k )} n [tj , t|] = {vr (x;)} (5.25) 
/or some i G {1, • • • , fc}, and 

{vr (xi) , • • • , vr (x fc )} n - o", ^ ] ^ 0, {vr (a*) , • • • , vr (x k )} D [t+,t+ + a) ^ 0, 
i/ien /or any optimal assignment map S G Map [£, fc], 

5- 1 (i) = {j:vr(y,)G [*2,4]}- 



Proof. For any j G S^ 1 (i), using ti = ^2 = tj m Corollary 5.2 and the fact vr (x j) G 
[tJ,oo), we have vr(yj) G [iJ,oo). Similarly, using t\ = tf, t2 = t£ in Corollary 5.2, and 
the fact vr(xj) G (— oo,tf], we have n (yj) G (— oo,t^]. Thus, 7r (yj) G ' ^i 1 ■ This shows 
that S-^i) C {j :vr( % ) G [tj,tf]}. 

On the other hand, for any j with vr (yj) G [i^ > n 1 > we have vr (yj) G [i^ > oo) and vr (yj) G 
(— oo, Using Corollary 5.2 again, we have vr (xsti)) G [i^ ,oo) and vr (a^sm) G (— 00, if]. 
Thus, vr (x su) ) G [ta ,4]. By g,S(j') = »■ Therefore, {j : vr (%•) G [tj,tf]} C S" 1 (i). 

□ 



6. State matrix 

In this section, we show that the properties of optimal assignment maps explored in 
previous sections can shed light on the search for those maps. The analysis is built upon a 
notion of state matrix defined as follows. 

Definition 6.1. Let U = (u s h) be an kx £ matrix with u s h G {0, 1}. The matrix U is called 

(1) a state matrix for an optimal assignment map S if S (h) ^ s whenever u s h = 0. 

(2) a uniform state matrix if U is a state matrix for any optimal assignment map. 

One could think of a state matrix as an information set of a planner during the search 
process for optimal assignment maps. An entry u s h = (or u s h = 1) simply denotes that 
the planner has (or has not) excluded the possibility of assigning household h to factory s. 
Recall that finding an optimal assignment map is to minimize the functional E a (S; x, b) 
over the set Map [£, k] whose cardinality is k . Any zero entry of a state matrix U for an 
optimal assignment map S may exclude as many as k^ 1 assignment maps in Map [£, k] 
from being S. The more zero entries in a state matrix U, the more information about S is 
contained in U. Consequently, we aim at finding a state matrix U for S with as many zero 
entries as possible, using properties of optimal assignment maps studied in previous sections. 
When U has exactly one non-zero entry in each column, S is completely determined by those 
non-zero entries in U. 

We first explore the implication of Theorem [2] on the search for optimal assignment maps 
in the context of state matrix. For any state matrix U = (uij), we consider a k x £ matrix 

W v = (wij (U)) 
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where 



W; 



(U) = p a (wi (U) , rij) with Wi (U) := ^ u ih n h , 



h=l 



and the function p a is given in (4.7). Here, Wj (U) denotes the maximum amount of com- 



modity produced at factory i one could conjecture using the existing information in state 
matrix U. 

For any state matrix U, define 

F s (U;nj) := \J{z£R m : \\z - Xl \\ < w sj (U) \\z - x s \\} 



and 



n s (U; nj ) :=p|{zGM m : \\z - x s \\ < Wij (U) \\z - Xi \\} , 



for any s = 1, • • • , k and j = 1, ■ ■ ■ ,£. By Lemma 4.1 each F s (U; rij) is the union of k — 1 
open balls 

\2 



_ x ■ K,(C/)) 2V ; ' 1 -(«;„■ (EO) 

and each set f2 s ([/ ; rij) is the intersection of k — 1 open balls 



F s (U;rij) = [J B \xi + V ~// m2 ~ Xg ) ' , , " , ,2 'I ~ 



n s (U; nj ) = {~] B lx s + - 



(^(u)) 2 . _ T] tgg (CO 



Example 6.1. The matrix 

U(°) = (uij) with Uij = 1 for any i and j 
is a uniform state matrix. Then, Wjj(o) = (wij (U^)) with 



Xj, Xa 



w. 



Pa (1,%) 



which is independent of i. Here, 

F s (V (0) ; rij) = F s (rij) and n s (u {0) ;rij) = tt s (rij) 



where F s (rij) and s (rij) are given in (4- 15) and (4- 16). 
Example 6.2. Let S 6 Map [£, k] be an optimal assignment map. Define 

Us = (u^) 



with 



1, ifS(j)=i 
0, else 



Then, W\j s = (w^j (Us)) with 



(6.1) 



(6.2) 



wij (U s ) = p a \ ^2 n h , rij J = p a (m (b») , rij) . 

K S(h)=i 
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F s (U s ; nj) = F s s (nj) and Q s {U s ; n 3 ) = Q s s ( nj ) 



where F s (rij) and Qg (nj) are given in (4-12) and (4-13). 

Definition 6.2. Given two k x i real matrices U = (uij) and U = (uij), we define 

(1) U > U if Uij > Uij for each i and j. 

(2) U>UifU>UbutU^U. 

Proposition 6.1. Let U and U be two state matrices for an optimal assignment map S. If 
U >U, then 

Wu < W D 

and 



x s en s (U;nj) <zn s [U;nj) and F s (U; nj ) ^F s [U;nj 



(6.3) 
(6.4) 



Proof. For each i, since U >U, 



IV; 



(U) = ^ U ih n h > ^2 ^ ihnh = W i( U 
h=l h=l 

Then, since p a (•,%•) is decreasing, it follows 

Wij (U) = p a (wi (U) , nj) < Pa[Wi[U) , nj 



W; 



U 



for each i and j. By definition, we have both (6.3) and (6.4). 



□ 



Let U be a state matrix for an optimal assignment map 5. By definitions of U<® in ffi.ll 

C/(°) >U>U S . (6.5) 



and Us in (6.2), it follows that 



Thus, by (6.4), we have 



m) C F s (U;nj)<ZF s (U s ;nj), 

n s (u {0) ;nj^ C fl s (U;nj) Qn s (U s ;nj) . 

These relations, together with Theorem [2j immediately imply the following proposition: 

Proposition 6.2. Let U = (u s j) be a state matrix for an optimal assignment map S G 
Map [I, ft]. For some s and j, 

(1) if Vj eF s {U-n 3 ), thenS(j)^s; 

(2) if Vj e Q s (U; nj), then S (j) = s. 

Corollary 6.1. Suppose U is a state matrix for an optimal assignment map S E Map [£, k]. 

Let f/W 



ii 



be a k x i matrix with 



0, ifyjCFHU; 



11 ; 



else 



(6.6) 



Then, U^ is also a state matrix for S with U > U^ . 



Proof. If Uy = 0, then either U{j = or j/j G F l (U; nj). In the first case, since U is a state 
matrix for S, by definition, S (j) ^ i. In the second case, by Proposition 6.2 S (j) ^ i. 
Thus, f/W is also a state matrix for S with U > U^\ □ 
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We now explore the implication of Theorem [3] on the search for optimal assignment 
maps. Suppose U is a state matrix for an optimal assignment map S. If u s h = for some 
h G {1, • " j ^} an d s G {1, • • • , k}, then for each j ^ h with nj < nh, we consider the set 

F s ' h (U; nj ) := {zei m : \\z -y h \\+A (U) < w sj (U) \\z - x s \\} , (6.7) 

where 

Pa (rih +rij,nj) 



A (U) = max 



Wij (U) 



\y h - Xi\\ : for i E {!,-■■ ,k} with u ih = 1 



Lemma 6.1. Let U and U be two state matrices for an optimal assignment map S. If 
U >U, then 



F s ' (U;nj) C F s ' h (U; nj 
for any s G {1, • • • , k}, h G {1, • • • , ^} u>i£/i u s /j = 0, and nj < for j ^ h. 



(6. 



Proof. For each i, if Uih = 1, then = 1 as £7 > U. By (6.3), we have Wij (U) < Wij (U 
Thus, 



A ( U 



= max < 

< max 

< max 



Pa {n h + rij,nj) 



( U 



Pa (n h + nj,nj) 
(U) 

Pa (n h + nj,rij) 
w^ (U) 



\\Vh - Xi\ 

\Vh ~ 
\Vh ~ Xi\\ 



for i with u. 



h 



1 



for i with Uih 
for i with Uj/j 



A (10- 



Consequently, (6.8) follows from (ra. 7to. 



□ 



As a result, by (6.5), 



F s ' h (U ; rij) C F s ' ft (£/; n,-) C F s,/l (U s ; n 3 ) = Ff (rij) , 



where F s ,h (nj) is given in (4.17). The following proposition and its associated corollary 
follow from Theorem O 

Proposition 6.3. Suppose U is a state matrix for an optimal assignment map S G Map [£, k] . 

If yj G F s,h (U; nj) for some h ^ j with u s h = and nj < nh, then S (j) ^ s. 

Corollary 6.2. Suppose U is a state matrix for an optimal assignment map S G Map [£, k]. 
Let = [u^j be a k x £ matrix with 

0, if yj G F s,h (U ; rij) for some h / j with u s h = and nj < ^ ^ 

U s j ^ 6 Is & 

Then, is also a state matrix for S with U > . 

We now explore the implication of Theorem [4] on the search for optimal assignment maps. 
Suppose U is a state matrix for an optimal assignment map S. For each i G {1, • • • , k}, let 

(U) = { Vj : mj = 1} U {xi} . 

Clearly, 

*i (U) 5 (U S ) = *<, 



■V 
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where is defined in (5.11). 
Now, for vr : R m - 



given in (5.1), we define 

R; 



(6.10) 



max{||z — p — vr (z) v\\ : z G \Pj (U)} . 
Without loss of generality, we may assume that 

®i(U) = {y jh :h = !,-■■ ,Ni}U{xi} 

with 

Trfoj-J <7r(y h ) < ■■■ <tt (y jN ^j . 

Proposition 6.4. Suppose U = (u s j) is a state matrix for an optimal assignment map 
S G Map [£, ft]. For each i G {1, • • • , k}, let h G {1, • • • , 2Vj} and i* G {1, • • • , ft}. If 

{ n (yjh+i) - n (%'J . 71 ( x i) - 71 (yjh)} > 2CR i + K ( x i*) - 71 (yjh) l> ( 6 - 12 ) 



(6.11) 



mm ■ 



where C and Ri are the constants given in (5.3) and (6.10) respectively, then S (jt) ^ i for 
any t < h. Similarly, if 

xmn{n(y jh ) - vr [yj h _ x ) ,^(yj h ) - n(xi)} > 2CR, t + |vr(x;«) - vr (y jh ) |, (6.13) 

then S (jt) 7^ i for any t > h. 



Proof. Assume (6.12) holds but S (jt*) = i for some t* < h. Without loss of generality, we 
may assume that 

71 (%>) = ma *U (Vjt) - t < h,S (j t ) =i} . 
Note that vr (yj t *) < vr (yj h ) < vr (x^ by (6.11) and (6.12). Then, by TheoremjiJ there exists 
G ^>i\ {yj t * } with vr (?/j t ») < vr (z) < vr (xi) such that 



a z 



< vr (z) - vr (y jt , ) < 2CR l + |vr (a*.) - vr (y jt , ) | . 



(6.14) 



By the maximality of vr (yj t *) and z G {yj t * }, we know vr (yj h ) < vr (z). Thus, by the 
ordering in (6.11), we have 



mm 



(6.15) 



Therefore, 



7T (z) - TT (y jtt ) 

vr (z) - vr + vr (y jh ) - vr ( %> ) 



> min {vr (y jh+1 ) - vr (y jh ) , vr (z;) - vr (y jh )} + vr (y jh ) - vr ) , by fl6.15| ) 

> + |tt (s f .) - vr (jyj | + vr ( % -J - vr ( %> ) , by fl6T2| ) 

> 2Ci2i + |7r(xi.)-7r(j/ it ,)|, 

a contradiction with (6.14). This proves (6.12). Similar arguments give (6.13). 



□ 



For each i = {1, • • • , k}, and A G K, denote 

Ii(A) := 

Then, for each vr : M m — > M given in (5.1 ), we define 



(— oo, A], if A < vr (xi) 
[A, oo), if A > vr (xi) 



F\ (U;nj) := {z G M m : vr (z) G I* (vr for some j h satisfying (6.12) or (6.13)} 
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Corollary 6.3. Suppose U is a state matrix for an optimal assignment map S £ Map [£, k\. 
Let = 



be a k x £ matrix with 



ii 



(3) 



if Uj £ F\ (U;nj) for some tt 
else 



(6.16) 



Then, is also a state matrix for S with U > . 



Remark 6.1. Depending on spatial locations of households and factories, for each fixed 
i £ {1, • • • , k}, the planner may choose tt to be one of the standard coordinate functions in 
W 11 , i.e. 7r (zi, • • • , z m ) = z t for some fixed 1 < t < m. In this case, (6.12) and (6.13) may 
be simply expressed in terms of coordinates of xi 's and yj 's. Another reasonable choice is 
to set tt (z) = (z — pi, Vi), where 



(Pi, Vi) £ argmin < max \\z 



p — (z — p,v) v\\ : p, v £ M m with \\v\ 



1 



This will minimize Ri given in (6.10), because the line passing through pi in direction Vi, 
i.e. 

{Pi +t Vl :t£R}, 
provides the least supremum norm approximation for (U) in M. m . 

Given a state matrix U for an optimal assignment map S, we have used results from 
previous sections to provide three updated state matrices j = 1, 2, 3, for U. The next 
proposition makes it possible to combine them together into a further updated state matrix. 

Proposition 6.5. Suppose U = (uij) and U = (v,ij) are two state matrices for an optimal 
assignment map S S Map[£, k\. Then, the matrix U = (uij) given by 

Uij = mm {u^ , Uij} for all i and j 

is also a state matrix for S. 



Proof. If Uij = 0, then either Uij = or Uij = 0. Both cases give S (j) ^ i. 



□ 



Proposition 6.5 says that one could deduce more information from any two existing state 
matrices regarding the optimal assignment map. Using this proposition, we immediately 
have the following corollary. 

Corollary 6.4. Suppose U is a state matrix for an optimal assignment map S £ Map [£, k]. 

For each i and j, define 

. f~(l) -(2) ^(3)1 
Uij = mm^yX/X/ V, 



, -(1) -(2) , ^(3) 

where u\- , u\- and u\^ are given in 
is also a state matrix for S with U >U 



6.6), 



6.9) and (6.16) respectively. Then, U 



6.4 



This idea of updating a state matrix U into another state matrix U as in Corollary 
can be implemented iteratively to obtain an even further updated state matrix. Given any 
initial state matrix U (e.g. U = as in (6.1)) for an optimal assignment map S. For 
each n = 0, 1, 2, • • ■ , define 

U n+ i = U n with Uq = U. 
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This gives a non- increasing sequence of k x £ matrices {U n } whose entries are either or 1. 
Hence, there exists an N > 1 such that 

U > Ui > ■ ■ ■ > U N -i = U N = U N+1 

We denote this Un as U* . Clearly, the matrix U* is still a state matrix for S with U > U* 
and U* = U*. 

This updated state matrix U* contains more information about S than the initial state 
matrix Uq because U* contains more zero entries. In some non-trivial cases as illustrated 
in the following example, U* may have exactly one non-zero entry in each column. In such 
a situation, U* completely determines the optimal assignment map S. 



Example 6.3. Let U be a uniform state matrix (e.g. U = as in (6.1)), and suppose 
that 

n\ > ri2 > • • • > ri£. 

If for each j = !,■■■ 



( 



U, S j — . 



\ 



F s (U; nj )U |J F s ' h (U;n,) 



V 



l<h<j-l 



(6.17) 



/ 



for some Sj G {1, 
assignment map. 



, k}, then S : {1, • • • , £} — > {1, • • • , k} given by S (j) = Sj is the optimal 



Proof. It is sufficient to show that for any optimal assignment map S, it holds that S (j) = Sj 
for any j. Indeed, for any s / Sj , if u s j = 0, then S (j) ^ s because U is a state matrix for S. 
If u s j = 1, then by assumption (6.17), either yj e F s (U; rij) or F s,h (U ; rij) for some h < j 

" S {j) / s. If Vj e F s ' h (U; nj ) for 
= Sh / s. In the later case, since 



6.2 



with Sh 7^ s. If yj G F s (U;rij), then by Proposition 
some h < j with Sh / s, then either S (h) / Sh or S{h 

yj € F s,h (U; rij) C Fg h (rij) and rij < rih, by Theorem pi we still have S (j) / s. Thus, in 
all cases for any s ^ Sj, we know 

either 5 (j) / s or S (h) ^ Sh for some h < j. (6.18) 

Consequently, when j = 1, we always have 5(1) ^ s for any s ^ s%, and thus S (1) = si. 
Using (6.18) again, we get S (2) ^ s for any s / S2, which yields S (2) = S2- Repeating this 
process leads to the conclusion that S (j) = Sj for any j £ {1, ■■■,£} ■ □ 



7. Conclusion 

This paper proposes an optimal allocation problem with ramified transport technology in 
a spatial economy. A planner needs to find an optimal allocation plan as well as an associated 
optimal allocation path to minimize overall cost of transporting commodity from factories 
to households. This problem differentiates itself from existing ramified transportation liter- 
ature in that the distribution of production among factories is not fixed but endogenously 
determined as in many allocation practices. It's shown that due to the transport economy 
of scale in ramified transportation, each optimal allocation plan corresponds equivalently 
to an optimal assignment map from households to factories. This optimal assignment map 
provides a natural partition of both households and allocation paths. We develop methods 
of marginal transportation analysis and projectional analysis to study properties of optimal 
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assignment maps. These properties are then related to the search for an optimal assignment 
map in the context of state matrix. 

The ramified optimal allocation problem studied in this paper provides a prototype for a 
class of problems arising in spatial resource allocations. One natural extension is to allow 
the locations of factories {x\,X2, ■ • ■ , Xk} to vary which then gives rise to an optimal location 
problem. An analogous optimal location problem in Monge-Kantorovich transportation has 
been extensively studied as in McAsey and Mou [IT] , Morgan and Bolton [20] and references 
therein. Meanwhile, one may consider another extension of the ramified allocation problem 
by generalizing the atomic measure b of households to an arbitrary probability measure \x, 
not necessarily atomic. In particular, when fi represents the Lebesgue measure on a domain, 
a partition of fi given by an optimal assignment map may analogously lead to a partition 
of the domain. This consequently gives rise to an optimal partition problem of dividing the 
given domain into k regions according to ramified optimal transportation. 
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